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Abstract 

An infinite system of nonlocal, individually confining solitons is considered 
as a model of high-density nuclear matter. The soliton-lattice problem is 
discussed in the Wigner-Seitz approximation. The cell size is varied to study 
the density dependence of physical quantities of interest. A transition to a 
system where quarks can migrate between solitons is found. We argue that 
this signals quark deconfinement. The model is applied to the calculation of 
selected in-medium properties. 
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I. INTRODUCTION 



High-density strongly-interacting matter is in the focus of attention of nuclear research 
for several reasons. Studying large (on the order of 100 fm^) volumes of dense hadronic 
matter experimentally is one aspect of the general effort to extend the investigation of the 
nuclear phase diagram beyond standard nuclear matter density and zero temperature. This 
experimental program is carried out at a set of accelerator facilities capable of colliding heavy 
nuclei at increasing energies. Such collisions provide the only way to access conditions in 
the laboratory that were dominant in an early stage of the evolution of the Universe and 
are relevant today in certain celestial objects and events, like dense stars and supernovae. 

It is of particularly great interest to identify and characterize in the laboratory the 
transition where strongly-interacting matter no longer appears as a collection of hadrons, 
but as deconfined quarks and gluons in an extended space-time domain. The high-density, 
low-temperature region of the nuclear phase diagram is especially important, as it may 
provide access to deconfined quarks without involving copious particle production and other 
effects of high excitation. Experimental results from the Brookhaven Alternating Gradient 
Synchrotron (AGS) indicate that this regime may be reached in the so-called full-stopping 
scenario achieved at the AGS. 

In the present paper we discuss a static approximation for this kind of high-density 
nuclear matter. Since we want to address the transition to the quark-gluon phase, we start 
with a description of the nucleon in terms of the underlying degrees of freedom. Ideally, such 
calculations should be carried out in the framework of quantum chromodynamics (QCD). 
However, as long as the solution of QCD at nuclear length and energy scales (nonperturbative 
regime) remains out of reach, modeling of QCD will play an important role. Here we focus on 
describing strongly-interacting matter at zero temperature, as a function of density. Soliton 
matter has been used to model high-density hadronic matter earlier. Bound states in 
a background soliton form the basic idea for the popular description of baryons with the 
Skyrme lagrangian. |Q 
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The elementary building blocks of the model proposed here are provided by the Global 
Color Model (GCM). The GCM admits soliton solutions with an intrinsically generated, 
extended qq meson field, in contrast to e.g. the Color Dielectric Model, |Q which uses an 
external field to generate a cavity (where quarks can propagate) in the vacuum. In addition, 
the individual GCM solitons are confining, as there are no poles in the quark propagator 
outside the region where the meson field is nonzero. The GCM observes the global 
symmetries of QCD, but is not locally gauge invariant. It has been used successfully to 
model low-energy QCD, as illustrated by the reproduction of chiral perturbation theory 
results, [pipip^ meson form factors, Jill and spectra, |jl2[ and both, the soliton (mentioned 
above), and the Faddeev ||13| description of the nucleon. A more exhaustive summary of 
these successes and of current work on hadron physics based on the GCM can be found in 
a recent review. [|14| 



To address high-density nuclear matter, we consider a lattice of GCM solitons. As this 
model maintains only global color symmetry, our work can be considered complementary to 
approaches that are concerned with a complete treatment of the color degrees of freedom 
for the description of nuclear matter. |TB[ Furthermore, the kinetic energy of the nucleons 



(solitons) is neglected in the present work. While this can be considered a reasonable 
approximation at the lowest temperatures, it implies that no quantitative agreement with 
e.g. the saturation density of nuclear matter should be expected. For the time being, 
we are more interested in identifying qualitative changes in the behavior of the system 
with increasing density, which can justify further work on extended strongly-interacting 
matter with this model. We describe the the soliton lattice at the mean field level, utilizing 



the Wigner-Seitz approximation. As a consequence, the present study is restricted to 
spherically symmetric mean fields. We are going to calculate excited states with higher 
orbital angular momenta in the spherical background field. 

The density of the system can be varied by changing the size of the Wigner-Seitz cell. This 
allows the study of stationary energies and in-medium properties as functions of the density 
in the model. Even though the simplicity of our picture precludes detailed quantitative 



predictions, we find a very interesting qualitative feature, tlie occurence of a transition from 
a color insulator to a color conductor at a certain density. We argue that this signals the 
deconfinement transition in the model. We show illustrative results, such as the axial- vector 
coupling constant and the correlation length for pion-like correlations as a function of density 
to highlight possible uses of the model. 

The paper is organized as follows. In Section |l| we review the hadronization of QCD in 
the context of the GCM to set the stage for the present work and to present the coupled 



equations defining our numerical problem. In Section |ITl| we discuss our choice of the gluon 
propagator and present results on a single soliton. Section |V| describes how we put the 
GCM "on the lattice," gives details on the solution method, and presents our results on 
soliton matter as a function of density, including energy bands and in-medium properties. 
We discuss the significance of these findings. Finally, in Section ^ we summarize the present 
status of the project and outline our future plans. 

II. FROM QCD TO THE GCM SOLITON EQUATIONS 

As quarks and gluons are directly unobservable, it is natural to seek a description of 
low-energy strong-interaction phenomena in terms of effective hadronic degrees of freedom. 
An example of the successes of such modeling is provided by Quantum Hadrodynamics. |T7 



Ideally, the effective degrees of freedom should be derived from the QCD Lagrangian. An 
approach to connect QCD and effective hadronic field theories can be formulated in terms 
of functional integral methods. The strategy is to transform the integration variables from 
quark and gluon fields to hadron fields. In order to make our discussion reasonably self- 
contained, here we review the major steps leading to the hadronic fields that play a central 
role in the present work. 

One particular implementation of the above ideas is in the framework of the Global 
Color Symmetry Model (GCM), which starts with a truncation of QCD, leading to the 
Euclidean action 
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= — J d'^x d'^y 



9^ 

g(x)(7 ■ d + m)6{x - y)q{y) + —f {x)Df,y{x - y)]t{y) 



Here = q{x)^'~^yq{x) is a local quark current, with Euclidean Dirac matrices 7,^ and 

Gell-Mann matrices A^. The two-point gluon function, -D^i,, can be considered the phe- 
nomenological input point for the model. Using a Feynman-like gauge, D^^ = 5^yD{x — y), 
the gluon propagator is particularly simple, and provides a single parameter function for the 
GCM. In the limit D{x — y)^ 5(x — y), the GCM reduces to the local Nambu-Jona-Lasinio 
model. [0 In ([I|), m is a current quark mass, which will be taken to be zero in the following. 
Our choice of the gluon propagator is dictated by the requirement of quark confinement, 
and will be detailed in Section The GCM has the global symmetries of QCD, but lacks 
local gauge invariance. 

The partition function in terms of the quark degrees of freedom can be written as 

Z = N jvqVq exp(^[gg]) , (2) 

where the functional integration Vq implies integration over all values of the quark fields, 
and is a normalization constant. To exhibit nonlocal quark- ant iquark structures in the 
action, a Fierz reordering may be performed, |TU| which transforms the current-current term 



as 



in® 

1 j d'x d'y j;{x)D{x - y)j;{y) = ~j d'x d'y J%x,y)D{x - y)j'{y,x) . (3) 

Here, J7'^(x, y) = q{x)A^q{y) can be looked upon as a bilocal current with a quark- ant iquark 
structure and quantum numbers specified by 9. The quantity is a direct product of Dirac, 
flavor, and color matrices, and contains both, color-singlet and color-octet terms. We focus 
on the color-singlet sector in this work, ignoring correlations that correspond to diquark 
degrees of freedom. 

To cast the partition function in terms of Bose fields, auxiliary nonlocal fields, B^{x, y), 
are introduced, and the partition function is multiplied by unity in the form 

l3%x,y)B%y,x 



1 = N' JvBexp - J d^xd'^y 



2g^D{x-y) 



(4) 



After the transformation B^{x, y) B^{x, y) + g^D[x — y)J^{y, x), the action is bihnear in 
terms of the quark fields and the Grassman integration can be performed. This yields the 
action in terms of bilocal Bose fields as 

S[B] = TvLn G-^[B] - f d'x d'y^-^^1^^^-^ , (5) 

J 2gW{x - y) 

where G'^{x, y) = (7 ■ <9 + m)5{x — y) + A^B^{x, y) is the inverse quark propagator. 

The replacement of the quark fields with Bose fields is, in principle, an exact functional 
change of variables. Observables calculated from the partition function are not affected by 
the variable transformation, but are now expressed in terms of the Bose degrees of freedom, 
provided the entire sum over 9 is kept. This is impossible in practice, and the truncation 
scheme used can be developed into a systematic method of approximation. To retain the 
chiral content of the QCD action, at least two Bose fields need to be kept (see below). 

The classical vacuum configuration Bq is identified by 5S/ 6B^ = 0. This produces a quark 
self-energy, T,{x — y) = A^Bq{x — y) satisfying a Schwinger-Dyson equation. In momentum 
space 

m = n ■ p|-4(/) - II + B(p^) = fj^,Di,- ^ ^ I ^ ^^^^ . (6) 

Numerical solutions for the amplitudes A{p^) and BijP') are now available at different 
levels of sophistication. pO[| As detailed in Section |ITT| , our choice for this explorative study is 



governed by simplicity within the context of the requirement of confinement. It is important 
to note that the amplitude B{p'^) plays a dual role in the model: it also acts as the distributed 
vertex for coupling the quarks to the qq Goldstone modes. [[7|,pT| 

The fluctuations &{x, y) = B^{x — y) — Bq are identified as the propagating Bose fields. 
If the color-singlet scalar-isoscalar and pseudoscalar-isovector fluctuations are retained, the 
formalism can be adopted to the requirement of chiral invariance by the variable transfor- 
mation 

A^^^(x,y) = ^x(i?)e^^-^(«)/^^, (7) 



where r = x — y and R = {x + y)/2, are relative and cm-like coordinates, respectively, f.„ is 
the pion decay constant, and it has been assumed that the on-shell form factor B can also 
be used off-shell. As a further simplification, the = point on the chiral circle can be 
fixed. In this case the radial fluctuations away from the chiral circle coincide with the scalar- 
isoscalar field variable prior to the transformation. In the numerical work that follows the 
single fiuctuation field x corresponding to this situation will be used, the notation serving as 
a reminder for possible genaralizations to restore chiral symmetry in the numerical model. 

Letting m — 0, the action (up to a constant) can be written as a sum of fermionic and 
bosonic terms: 

= -Tr[ lnG-i(/i,x) -lnG-i(0,x)] + / d'R [^(^^x)' + f/(x')] • (8) 

The chemical potential (/i) dependence of the fermion term in equation (H) ensures that a 
meson source from the valence quarks will be generated. |^ The U{x^) term is the effective 
meson self-interaction. For n = the inverse quark propagator takes the form 

G-\x, y) = ^-d, A{x -y) + f-^B{x - y)x (^) , (9) 



and the saddle-point configuration turns out to be x = /^r. 

Since y) is time-translationally invariant, stationary eigenstates of the form Uj{x) 

can be obtained from a self-consistent Dirac equation, which in momentum space takes the 
form 



[.7 ■ pA{p') + i?(/)]«,(p) + f-' J j^,B 



q\2 



x(p-qK(q)=0, (10) 



where p and q are the three-momenta corresponding to the center-of-mass and relative 
variables, respectively. Note that X = X ~ fir, and as ^4 = ^4 = iej, where ej is the energy 
eigenvalue, the meson vertex B has an energy dependence. It can also be seen that a 
wave-function renormalization appears with the renormalization constant Zj given by ^ 



The meson field equation = may be summarized as 



Vx(z) + 



6U 



+ Qx(z) = 



with the meson source provided by the valence quarks according to 

"x + y 



(12) 



(13) 



Equations (|TOD and (|T2D form a system of coupled differential equations for the quark wave 
functions and the meson field, which need to be solved selfconsistently, with the appropriate 
boundary conditions. The different boundary conditions distinguish the single-soliton case 
from a lattice of solitons. 



III. SINGLE SOLITON 

A. Gluon Propagator 

One can take the point of view that the gluon propagator, D^^ in Eq. (|l|), represents 
all phenomenological input to the model. With the appropriate choice of D, the GCM can 
reproduce key features of QCD, such as confinement. Phenomenologically successful early 
work with the GCM ^ employed a delta-function gluon propagator, 

g^D{q) = 37iW6'^^\q) . (14) 

The momentum-space delta function turns the rainbow approximation Schwinger-Dyson 
equation into an algebraic equation, and the Munczek-Nemirovski quark propagator p3|] 
results: 

A{p') = { _ , B{p-) = {' \ ■ (15) 

Its simplicity, in addition to its confining nature, make this form of the gluon propagator 
particularly appealing. It has a single strength parameter, a, which (for a fixed energy) 
controls the spatial extension of both A and B: the larger a, the more localized A and B 
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are in coordinate space. The lack of solutions to the equation + M'^ijP') = 0, (where 



M = B/A is the dynamic quark mass) indicates that (|T5|) produces quark confinement; 
there is no on-mass-shell point, thus the propagation of a quark in the normal vacuum is 
prohibited. ^ 

Though phenomenologically successful, the above simple quark propagator encounters 
difficulties when improvements are attempted regarding analyticity. Partly for this rea- 
son, recent studies in the hadronic sector moved away from the point of view of providing 
the input at the level of the gluon propagator and, accordingly, from the simple form (|14|). 
In order to focus on nuclear matter and to keep complications from the GCM to a minimum, 
we use the delta-function gluon propagator (|^) throughout this work. Note also the recent 
observation [jl4[ that the parametrization (|l^) may well express the infrared behavior of the 



gluon two-point function in QCD. [25 



B. Single Soliton Results 



With the self-energy amplitudes A and B of (|T5D, the coupled equations (0) and (|T2D 
can be solved self-consistently. For convenience, the Dirac equation is solved as a matrix 
equation in momentum space, using Gauss- Legendre quadrature, while the Klein-Gordon 



equation is solved in coordinate space using a functional version of Newton's method. |22 



In order to check our numerical procedure, we have satisfactorily compared our results on 
the three valence-quark soliton to the results of an earlier investigation. [0,^] In addition, 
we carried out calculations for a range of the parameter a. 

Assuming spherical symmetry, the quark spinors can be decomposed as 



■u(r) 



(16) 



g{r) 
ia ■ kf{r) 

Here, y'-l^ (f) is a vector spherical harmonic, f represents a unit vector in the direction of r, 
and we have suppressed the quantum-number labels on /, g, and u. The numerical task is 
now reduced to the calculation of the radial functions g{r) and f{r). 
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The radial parts of the upper and lower components of the quark wave function and the 
meson field (x) for a = 1.04 GeV are shown in Fig. 1. These results should be compared 



to the results of where a = 1.04 GeV was chosen to fit the experimental value of 
/vr. For our nuclear matter studies, we find it more important to have a reasonably close 
correspondence to the root mean square charge radius of the proton. As discussed in Ref . , 
this does not fix a in lack of explicitly considering the pion field, but it appears to call for 
larger values of the strength parameter. With a = 1.35 GeV we get (r^)^/^ = 0.67 fm for 
the RMS charge radius of the proton. Assuming that the pion cloud will increase this value 
by about 25-30 %, we get close to the experimental value. To leave room for the uncertainty 
of this estimate, we consider the range 1.04 < a < 1.45. In the following figures, when a is 
fixed, we display results for a = 1.35 GeV. In Fig. 2 we show the Dirac wave function and 
X for a single soliton with this value of the strength parameter. 

Figures 1 and 2 represent typical results for a single soliton. It is noteworthy that the 
large and small components of the quark wave functions decay faster than exponential, ||2^ 
and are essentially zero by r = 3 fm for the values of a considered. The range of x decreases 
and the magnitude of x at r = increases as a increases. The changes of the x fi^ld can be 
attributed to the decrease in the range of the distributed quark-meson vertex B in coordinate 
space as a increases. The narrower B, the closer x and y have to be in Eq. (|13D to give a 
larger value of the source term in ([l2|) , which in turn infiuences the gradient of the x field: a 
larger source term leads to a stronger gradient (assuming no change in the self-interaction). 
In Fig. 3 we plot the dependence of the quark eigen energy e on the single input parameter 
a. The increase of e can be associated with the increasing absolute value of x at r = 0. The 
spreading of x in momentum space causes the coupling term in (|10|) to decrease, making the 
valence quarks to be less tightly bound. 



IV. SOLITON MATTER 
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A. Wigner-Seitz approximation 



As a means of describing nuclear matter, we consider an infinite collection of solitons. At 
the lowest energies the solitons are expected to arrange themselves in a crystal lattice. ||2^ 
Accordingly, the single-quark eigenenergies will develop into energy bands. For simplicity, 
we assume a simple cubic crystal (sc). For a periodic lattice, the Dirac wave function must 
be invariant to a lattice translation, so the solutions must have the Bloch form 

<(r)=«^(r)e^'^■^ (17) 

where m is the lattice momentum and Mm(r) is a Dirac spinor which has the periodicity of 
the lattice. To solve for the Bloch functions we employ the Wigner-Seitz approximation, fl^ 



This amounts to considering a spherical cell of radius R and solving for m = in ([T7|). The 
full solution Mm (r) is then approximated by no(r) e*™ ' ^ . Changing the density will be 
implemented by varying the cell radius R. 

The Wigner-Seitz approximation places boundary conditions on the Dirac spinors. These 
conditions express the requirement that the upper component of the wave function must be 
periodic and anti-periodic for the bottom and the top of the band, respectively. We focus 
attention on the lowest-energy state of the band, for which the above, together with the 
r = boundary conditions, imphes 

^?V)Uh = /WU = o, (18) 

where g{r) and /(r) represent the radial parts of the upper and lower components of the 
Dirac wave function (|16]), respectively. In addition, the meson field solution of eq. (|12|), 



which also appears in the source term of the Dirac equation (|ToD must now be periodic in 
r, so that 

x{r + 2R)=xir); x'{r)l=R = 0, (19) 
where x(r) is the radial part of the meson field. 
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B. GCM on the Lattice 



It is convenient to solve the Dirac equation in momentum space, while the nonlinear 
Klein-Gordon equation is easier to handle in coordinate space. The Dirac equation (p!OD can 
be written in coordinate space as 







1 



d'y (-746,- + 7 ■ V)A(x - y) + 5(x - y) + — fi(x - y)x 



x-y 



u,{y). (20) 



We seek solutions with the boundary conditions ([T^)- These can be incorporated using a 
Fourier expansion of the form 



oo oo oo 



/W= E E E nK 

n\=—oo 712=— oo 713=— oo 



,ik„ ■ X 



(21) 



where k„ = ^ is a wave number vector, with n = {ni, ?t,2, "^,3}. Expanding A, B, u, and the 
meson-field source, we integrate over y and use orthonormality to get an equation for the 
Fourier components of the j*^ quark wave function 



-eA(^kn) + B(^kn))Snm ~knA{kn)Snm 

knA{kn)Snm (c^(^n) + B{k,„))5n 



+ 



47r 



X]m,=0 ^mYoi^n^ k„ 







J2m=0 ^m^li^ri} km) 

where gn = g{kn), fn = fiK), and 

k -I- k ' 



9m 
fm 



(22) 



Vi 



B 



x(k„ - km)P;(cos6')d(cosi 



(23) 



To get this final form we have written the Dirac wave functions as 



■ufk^ 



(24) 



g{.kn) 

ia ■ KfiK) 

and used the spherical symmetry of the Wigner-Seitz cell. If m = 1, 2, . . . , M , then equation 
(E^) is an 2M-by-2M eigenvalue problem for the energy eigenvalue e. The quantity B/A 



12 



plays the role of a dynamic mass and the scalar part of the self-energy B also acts to couple 
the quarks to the meson field via (0). The self-energy terms have an e dependence which 
makes this a highly non-linear problem. The Vi term represents the Legendre coefficient for 
the meson field in the presence of the distributed coupling B. One needs to solve the Dirac 
equation and the Klein- Gordon equation for the meson field (0) self consistently. 



C. Details of Solution 

To solve for a soliton lattice, we first pick a starting meson field and search for the lowest 
energy eigenvalue of Eq. (]22|). Technically this means finding the energy e which makes the 
determinant of equation (p^ vanish. We start at e = and work upwards in energy until the 
determinant changes sign. We then use the bisection method to find the root. Care must be 
taken so that the initial steps are sufficiently fine in e not to miss the lowest root. With the 
root in hand, we can solve for the Fourier components of the Dirac wave functions. For this 



we first perform a lower upper triangular (LU) decomposition and use inverse iteration. p8 
The momentum-space meson-field source term is constructed from the Dirac wave functions 
and we transform the source to coordinate space for use in the nonlinear Klein-Gordon 
equation for the meson field (|12|) . To solve this nonlinear equation, we treat equation (p!2[) 
as a functional of x and use Newton's method. Once the Klein-Gordon equation is solved 
for the new meson field, we start over with the Dirac equation in this modified meson field. 
We iterate until convergence of the quark wave functions is achieved, which takes between 
three to six iterations to reasonable accuracy. 



D. Lattice Results for the Fields 

When the convergence of the Dirac and Klein-Gordon equations has been reached, we 
have the selfconsistent ground state quark wave functions, meson field, and quark energy at 
our disposal. For i? = 10 fm and a = 1.35 GeV a quark energy e = 530 MeV is obtained. 
This reproduces the single soliton value (e = 537 MeV) within 2%. Fig. 4 displays the upper 
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and lower components of the Dirac wave functions along with the meson field for i? = 5 fm 
(with a = 1.35 GeV). This still is a relatively large cell size and the results resemble the 
single-soliton case in shape, but the wave functions and the meson field go to zero faster 
than for a single soliton (compare to Fig. 2). The wave functions are pushed inwards by the 
boundary: each soliton is isolated around the center of the cell. The quark eigen energy for 
this case is e = 526 MeV, slightly smaller than at i? = 10 fm. 

For Fig. 5, we have decreased the size of the Wigner-Seitz cell to R = 1.5 fm while 
keeping the strength parameter a unchanged. Now the value of the meson field is different 
from zero at the edge of the cell. The upper component of the Dirac wave function, g{r), can 
no longer decay to zero within the cell and obey the boundary condition ([18|). Its value is 
substantial at the cell boundary. The lowest quark energy begins to increase to e = 542 MeV. 

Fig. 6 shows a cell with radius R = 0.85 fm and strength parameter a = 1.35 GeV. 
The meson field is non-zero at the cell boundary and g{r) is becoming relatively flat. As 
the quarks now strongly feel the presence of the neighboring cells, they become less tightly 
bound, and the eigen energy increases further to e = 621 MeV. Comparing Figures 4, 5, 
and 6, we conclude that there is a systematic evolution of the solutions as the cell size gets 
smaller (i.e. the density becomes larger) for a fixed value of the strength parameter. One 
important feature is that the quark distribution near the cell boundary becomes larger with 
increasing density: the solitons in neighboring cells begin to communicate. To clearly display 
this trend, we introduce the dimensionless variable r/R, and plot the upper component of 
the quark wave function and the meson field near the cell boundary as a function of this 
quantity in Fig. 7. This normalized variable is best suited for comparison between different 
values of R. The relative increase of the size of the meson field and of the large component 
of the wave function at the edge of the cell is now obvious. 
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E. Energy Bands 



Each soliton of the lattice contributes one level to each energy band. In the Wigner- 
Seitz approximation we need to calculate only the energy for the case when the crystal 
momentum m in equation (|T^ is zero. To approximate the top of the energy band, one may 
solve the equations with antiperiodic boundary conditions, P^JSIH ] or use an estimate for 



the band width. At the level of our computational accuracy we do not expect these two 



methods to give significantly different results and follow the simpler band-width calculation. 
We approximate the Dirac-Bloch wave function for an arbitrary crystal momentum as 

<*(r)=«o(r)e^^■^ (25) 

and use (^Sf ) to calculate the expectation value for the square of the Dirac Hamiltonian to 
estimate the lattice momentum dependence of the energy levels in the band as 



~bot 



+ m2]2, (26) 



where ebot is the energy of the bottom band. 

To obtain the possible values of m the lattice structure needs to be specified. For a 
simple cubic crystal of solitons and sides of length L = 2RN, the allowed values of the 
component of the lattice momentum in the direction of any of the three axes are 

271 , 47r Ntt 
m = 0, ±— , ±^---,— , (27) 

with the top of the energy band corresponding to m = ^ = Thus for the top energy 
band we obtain 

eto,= [eL + {^r]- . (28) 
We have performed the same estimate assuming body-centered and face-centered cubic 



lattices, the face-centered cubic deviating the most from the estimate p8|). Different as- 
sumptions in the lattice structure introduce an uncertainty of roughly 8 % in our results for 
the top of the band. 
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In Fig. 8 we show the three lowest energy bands for a = 1.35 MeV as a function of the 
density. The different symbols represent the calculated energies of the bottom of each band. 
On one point we indicate a typical uncertainty we associate with our computation. The 
main source of this uncertainty is the freedom in prescribed tolerances at different stages 
of the calculation. The lines across the symbols represent a polynomial fit to guide the eye 
and to facilitate the calculation of the top of the bands with the approximation (^Sj) . 

The lowest band (/ = 0, /' = 1 and j = |) is labeled lsi/2. The next lowest band has 
nonzero orbital angular momentum in the large component of the Dirac wave function (/ = 1, 
/' = 2, and j = |) and is labeled lp3/2- The next band is again an s-state, corresponding to 
a radial excitation. For very low density {R — > oo) the energy bands shrink to single levels 
and in the limit reproduce the energies of a single soliton (discussed earlier). As the density 
increases the ground-state band develops a minium. The low-density attraction between the 
solitons is a consequence of the boundary conditions on the quark wave functions (|18]). In 
particular, the upper component of the quark wave function is forced to have less curvature 
than in the case of a single soliton, leading to a lower value of the quark kinetic energy. 
At higher densities, where the solitons and the quark wave functions begin to overlap, 
the resulting repulsion overcomes the attraction and the ground-state quark energy starts 
to increase. In a more complete calculation one would like to attempt to fit the energy 
minimum to the saturation density of nuclear matter. This, however, requires to go beyond 
the present mean field treatment which lacks the details of the nucleon-nucleon interaction 
and nucleon kinetic energies. 

As the density increases, the top of the ground-state band approaches the bottom of the 
next unfilled band, and at ~ 2.6 po the highest energy state of the occupied ground state 
band intersects the bottom of the empty lp3/2 band. At this point it becomes energetically 
favorable for the quark in the highest-energy state to move into the empty "conduction" 
band. The system goes through a transition very similar to the insulator-conductor tran- 
sition in metals, and color conductivity sets in. Since in the new phase quarks are free to 
migrate from soliton to soliton, we identify this transition with quark deconfinement. ^ 
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Note in this context that we use uniformly filled bands in our calculation. Partial filling of 
the lowest energy band will increase the critical density in the model. 



F. In-medium properties 

The surrounding dense nuclear matter may significantly change single-particle properties 
like masses and widths, as well as coupling constants, cross sections, and other features 
relevant for transport modeling. For example, dilepton production experiments at CERN 
| ?T| seem to indicate a shifting p-meson mass in medium. Here we calculate the axial 



vector coupling constant in the model to illustrate how the density dependence of physical 
properties can be obtained. As another example of an in-medium property, we present the 
calculation of a correlation length. 

The axial vector coupling constant can be calculated in momentum space. Since there 
are no explicit pion fields in the model, only the valence quarks contribute. In a simple 
approximation 



9A = iydp p'A{-e' + p') [g\p) - ^f{p)^ , (29) 

where the | factor is due to the summation over spin and fiavor degrees of freedom. In 
Fig. 9 we plot qa as a function of density for several values of the strength parameter a. 
[The line serves to guide the eye.] For large cell size (low density) the single soliton result 
is approached. In general, the value of qa decreases with increasing density. This can 
be understood in terms of g{p) becoming more and more localized in momentum space as 
g{r) spreads out in coordinate space with decreasing cell size. The decrease continues until 
the critical density, which is the highest density up to which we trust our calculations based 
on a single cell. At the transition we physically expect a sudden increase in the value of qa 
to higher than its free-space value, as a consequence of the fact that the deconfined quarks 
sample a larger volume of phase space. The trend at higher densities that qa is smaller 
for smaller a is consistent with the narrowing of A in momentum space as the strength 
parameter decreases. 
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As another illustration of the utility of the model, we construct the pion-like bilocal 



fluctuation field q^y^'y^q, and calculate the current-current correlation function |^| 
p/ ^ ^ _ ( ^(ri)7/.75^(r2) u{v2)lf,l-,u{ri) ) 

To reduce the number of variables we set ri = r and r2 = 0. Averaging over angles and 
using the properties of the Dirac spinors makes it possible to write the correlation function 
as a function of one variable, the relative distance r. The result is 

r(r)= f . , (31) 

To extract a correlation length Rcor, (|3T|) is transformed to momentum space and we follow 
the methods usually applied in Hanbury-Brown and Twiss (HBT) types of analyses for 
bosons: the width of the momentum-space peak is inversely proportional to Rcor- 

In Fig. 10 we show the correlation length Rcor as a function of density for three values 
of a. [The line serves to guide the eye.] With increasing density the cells become smaller, 
and Rcor decreases. At the transition density a sudden increase of the correlation length is 
expected, as the quarks become free to migrate from cell to cell. For small density (large 
cell size) there is a large spread in Rcor as a function of a. This can be connected to the 
increasing spatial localization of the solitons as a increases. For smaller values of a the 
solitons spread out more, so that the lower component /(r) peaks at a larger distance away 
from the origin. For small cell size (large density) the wave function has no room in the cell 
to spread out, and Rcor becomes more independent of a. 

The above examples serve to illustrate how the model can be used to discuss the density 
dependence of in-medium properties. We plan to calculate other hadronic observables, like 
in-medium masses and further correlation functions in the future. 



V. SUMMARY 



We presented a generalization of the Global Color Model (GCM) to the many-soliton 
situation. In particular, the density dependence of the properties of an infinite system of 
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nonlocal, confining solitons was studied in the Wigner-Seitz approximation. We found that, 
at a critical density, an infinite system of solitons exhibits a transition from a geometry with 
one soliton localized at the center of each cell to a more uniform situation where quarks can 
migrate across cell boundaries. We argued that this transition signals quark deconfinement 
in the model. We have also calculated the density dependence of the axial vector coupling 
constant and of a correlation length as examples of in-medium properties. 

It should be kept in mind that the least elaborate one-parameter version of the GCM was 
used throughout this work. More realistic parametrizations of the quark self-energy functions 
should improve the accuracy of the description. In particular, the pion decay constant and 
the root mean square charge radius of the proton could be fitted simultaneously with a 
couple of parameters instead of just one. Such improvements in the hadronic sector can 
be incorporated in the model if a closer correspondence to experimental data is desired. 
Developments along these lines require explicit pion degrees of freedom in the model for 
inclusion into the calculation of the root mean square proton radius. The inclusion of 
explicit pions will also lead to the restauration of chiral symmetry in the model, as discussed 
following Eq. ([^). As this would open the way for chiral calculations, it promises to be an 
interesting line of future development. 

It would also be of interest to perform direct comparisons to QCD calculations. This 
can shed more light on the nature of the transition we found in the model. We consider 
the existence of this transition in the GCM-based soliton lattice to be our most important 
finding so far. We believe that this feature is sufficiently robust to survive in more refined 
versions of the model. 
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FIGURE CAPTIONS 



Fig. 1: Upper and lower components of the quark wave function for a single soliton with 
the corresponding meson field plotted for a = 1.04 GeV. The quark eigen energy is 
e = 356 MeV. 

Fig. 2: Upper and lower components of the quark wave function for a single soliton with 
the corresponding meson field plotted for a — 1.35 GeV. The quark eigen energy is 
e = 537 MeV. 

Fig. 3: The dependence of the quark eigen energy e on the strength parameter a for a single 
soliton. 

Fig. 4: Upper and lower components of the quark wave function for the sohton lattice with 
the corresponding meson field for R = 5.0 fm, a = 1.35 GeV. The quark eigen energy 
is e = 526 MeV. 

Fig. 5: Upper and lower components of the quark wave function for the soliton lattice with 
the corresponding meson field for R — 1.5 fm, a — 1.35 GeV. The quark eigen energy 
is e = 542 MeV. 

Fig. 6: Upper and lower components of the quark wave function for the soliton lattice with 
the corresponding meson field plotted for R = 0.85 fm, a — 1.35 GeV. The quark 
eigen energy is e = 621 MeV. 

Fig. 7: Upper components of the quark wave functions and the corresponding meson field 
for the lattice case plotted against the dimensionless variable r/R for R — 1.2, 1.3, 
and 1.5 fm, and a = 1.35 GeV. 

Fig. 8: The bottom and top energies for the lowest energy bands of the soliton lattice as a 
function of the density in terms of standard nuclear matter density (po = 0.17 fm~^), 
for a = 1.35 GeV. The symbols represent the calculated energies of the bottom of the 



22 



three lowest energy bands. For the top of the energy band we use (|2^) as an approxi- 
mation. An illustrative error bar is included. The curves represent a polynomial fit to 
guide the eye. 

Fig. 9: The axial vector coupling constant as a function of the density for a = 1.25 GeV, 
1.35 GeV, and 1.45 GeV. 

Fig. 10: The correlation length for pion-like currents as a function of the density for a = 
1.25 GeV, 1.35 GeV, and 1.45 GeV. 
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